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SCALING VIOLATIONS OF QUARK AND GLUON JET 
FRAGMENTATION FUNCTIONS IN E+E ANNIHILATIONS 
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Flavour inclusive, udsc and b fragmentation functions in biased and unbiased jets 
and gluon fragmentation functions in biased jets are measured in e+e" annihila- 
tions using OPAL data collected at energies y/s = 91.2 and 183-209 GeV. 



1 Introduction 



> 

lO ' Hadron production in high energy coUisions can be described by parton showers, 

\l foUowed by the formation of hadrons which cannot be described perturbatively. 

Gluon emission, the dominant process in parton showers, is proportional to the 
colour factor associated with the coupling of the emitted gluon to the emitter. 
. These colour factors are Ca — 3 when the emitter is a gluon and Cp =4/3 when 

' it is a quark. Consequently, the multiplicity of soft gluons from a gluon source is 

. (asymptotically) 9/4 times higher than from a quark source pp. The inequality 

' between Ca and Cp plays a key role in the explanation of the observed differences 



between quark and gluon jets: compared to quark jets, gluon jets are observed 
Qj ' to have larger widths [21, higher multiplicities |2I8| . softer fragmentation functions 

. |2l4l5j , and stronger scaling violations of the fragmentation functions [Sj . 

The fragmentation function, D^{x, Q^), is defined as the probability that parton 
a, which is produced at short distance, of order l/Q, fragments into hadron, h, 
' carrying the fraction x of the momentum of a. In this study, the momentum 

^ , fraction is defined as — Eh/Ej^t, where Ef^ is the energy of the hadron h and 

Ejet is the energy of the jet to which it is assigned. 

Three-jet qgg events are selected by applying a jet finder. Different jet finders 
result in different assignments of particles to jets: thus jets defined using a jet find- 
ing algorithm are called biased. In contrast, quark and gluon jets used in theoretical 
calculations are usually defined as inclusive hemispheres of back-to-back qq and gg 
final states, respectively. The hemisphere definition yields a so-called unbiased jet 
because the jet properties do not depend on the choice of a jet finder. Measure- 
ments of unbiased quark jets have been performed at many scales [61718) . Direct 
measurements of unbiased gluon jets are so far available only from the CLEG "W and 
GPAL |4I10| experiments, however. Recently, the GPAL experiment has measured 
properties of unbiased gluon jets indirectly |llll2j . 

In our study, we present measurements of quark and gluon jet fragmentation 
functions at = 91.2 GeV and ^/s = 183-209 GeV. The data were collected 
with the OPAL detector at the LEP e+e" collider at CERN. We measured seven 
types of fragmentation functions: those from biased as well as unbiased flavour 
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inclusive, udsc and b jets, and from biased gluon jets. While the two types of 
flavour inclusive jets have been measured many times, data on the other types of 
fragmentation functions arc still rather scarce. 

2 Analysis procedure 

The selection of hadronic Z and Z* /j* events and any other details of this analysis 
are described in 13 . In the inclusive hadronic event samples, we use the unbiased 
jet definition where the jets are defined by particles in hemispheres of the qg system. 
In the three-jet samples, we apply a jet algorithm and thus work with biased jets. 
Three jet algorithms are used: the Durham I 1 . Cambridge [13! and cone [T^ 
algorithms. The jet algorithm is forced to resolve three jets per event. The jet 
energies and momenta are then recalculated by imposing overall energy-momentum 
conservation with planar massless kinematics, using the jet directions found by the 
jet algorithm. 

The measured fragmentation function is defined here as the total number of 
charged particles, Np, in bins of xe and scale Q normalized to the number of jets, 
NjetiQ), in the bin of Q: 

1 dNp{xE,Q) , . 

NjetiQ) dxE ^ ' 

To measure the scale dependence, it is necessary to specify a scale relevant to 
the process under study. For inclusive hadronic events, the scale is y/s. For jets in 
three-jet events, neither y/s nor Ejet is considered to be an appropriate choice of 
the scale QCD coherence suggests ^H] that the event topology should also be 
taken into account. Similarly to previous studies, the transverse momentum-like 
scale, Qjet = -Ejet sin(i?/2), is used where is the angle between the jet with Ejet 
and the closest other jet. 

In this analysis, three methods are used to identify quark and gluon jets: the b- 
tag and the energy-ordering methods in biased three-jet events, and the hemisphere 
method in unbiased inclusive hadronic events. In addition, b tagging is used to 
separate udsc and b quark jets from each other, both for the biased and unbiased 
jet samples. The energy-ordering method only allows flavour inclusive quark jets 
to be distinguished from gluon jets. 

All three methods are in detail described in ^S]- In three-jet events, any of 
the three jets is used to extract the fragmentation functions. The purity of the 
b-tag (gluon) jet sample in LEPl data is estimated to be 90% (84%) and the 
efficiencies are 23% (40%). The corresponding purities in LEP2 data are 60% 
(80%) and efficiencies 27% (45%). The purity of the b-tag hemispheres in LEPl 
data is estimated to be 99.7%, while in LEP2 data it is 75%. 

After subtracting the remaining background from LEP2 data using MC es- 
timates, the data and MC distributions at detector level are multiplied by corre- 
sponding inverse matrices to get the distributions at level of pure quarks and gluons. 
As a last step, the data are corrected for effects of limited detector acceptance and 
resolution using correction factors obtained from MC events. 
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3 Results 

The measured fragmentation functions are presented with emphasis on the scale 
dependence (the xe dependence and other comparisons are shown in and are 
compared with previous measurements as weU as with theoretical NLO predictions. 
The term scale in the following figures stands for Qjet in case of biased jets and 
^/s/2 in case of unbiased jets. The published unbiased jet results and the NLO 
predictions are scaled by ^ since they refer to the entire event, thus to two jets. 

In Figs. la)-d) the results for the udsc, b, gluon and flavour inclusive jet frag- 
mentation functions are presented. The LEPl unbiased jet data correspond to 
^/s = 91.2 GeV. The LEP2 unbiased b jets are measured in y/s range of 183- 
209 GeV, while the unbiased udsc and flavour inclusive jets are measured in three 
y/s intervals: 183-189, 192-202 and 204-209 GeV. The quark biased jet data from 
LEPl cover the region Qjet = 4-42 GeV, while those from LEP2 cover the region 
Qjet = 30-105 GeV. The results from the region 0.01 < xe < 0.03 are not shown 
but they are presented in J^l- As also discussed in [121, the results are found to 
be consistent with previous measurements. The udsc jet results above 45.6 GeV, 
the gluon jet results above 30 GeV (except for the gind jets, see and the b jet 
results at all scales except 45.6 GeV represent the first measurements. 

The data are compared to three theoretical predictions: Kniehl, Kramer and 
Potter (KKP) [Tn|, Kretzer (Kr) gOI and Bourhis, Fontannaz, Guillet and Werlen 
(BFGW) [21]. For the udsc jet fragmentation function (Fig. la)), all three predic- 
tions give a good description in the entire measured phase space, except for the 
lowest Xe bin where the KKP calculations overestimate the data, and the highest 
Xe bin where the data are underestimated by the Kr and BFGW calculations. 

The situation is rather different for the b and gluon jet fragmentation functions 
(Figs, lb) and Ic)) where the description of the data by the NLO predictions is 
worse and where there are significant differences between individual NLO results, 
the latter being expected due to differences in the calculations as discussed in [O] . 
In Fig. lb) the KKP prediction is deficient with respect to the data for xe > 0.12. 
In this context, it is important to note that B hadron decay products are indirectly 
included in theory predictions since they are present in the data sets to which the 
fits were made. 

For the gluon jet fragmentation functions, the two alternative methods of iden- 
tifying gluon jets described above are examined, see Fig. Ic). The Qjet binning is 
not the same for the two methods because of their different regions of applicability. 
A satisfactory correspondence between the b-tag and energy-ordering methods is 
found in the entire scale range accessible. The data tend to show larger scaling 
violations than predicted by any of the calculations. 

The results for the flavour inclusive jet fragmentation functions, presented in 
Fig. Id), are seen to be consistent with published unbiased jet data from lower 
energy e+e" experiments (TASSO, MARK II, TPC and AMY) and previous 
OPAL results 22 25 24J. The data are also compared to the NLO predictions of 
KKP, Kr and BFGW which all give a reasonable description of the data in the 
region of 0.06 < xe < 0.60 and over the entire scale range. 

A good correspondence found between the results from biased and unbiased jets 
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in all four figures suggests that Qjet is an appropriate choice of scale in three-jet 
events with a general topology. A similar conclusion was previously made in 
The MC study presented in ^^1) however, demonstrates that the bias introduced in 
the gluon jet identification is not negligible for xe > 0.6. In each of these figures, 
the scaling violation seen in the data is positive for low xe and negative for high 
xe- It is more pronounced in the gluon jets than in the quark jets. 
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Figure 1: Scale dependence of fragmentation functions in different x-^ bins. The scale denotes 
Qjet for the biased jets and y/s/2 for the unbiased jets. The inner error bars indicate the statis- 
tical uncertainties, the total error bars show the statistical and systematic uncertainties added in 
quadrature. The data are compared to the NLO predictions by KKP Q^, Kr and BFGW 1^ . 
a) udsc jets, b) b jets, c) gluon jets from b-tag (BT) method and energy-ordering (EC) method. 
In addition, the published gind j^ts |1] are shown at E^^t = 40.1 GeV. d) flavour inclusive jets. 
In addition, the published unbiased jet data by TASSO, TPC, MARK II. AMY W and OPAL 
mVimH are shown. 



